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The dynamics of the population imbalance of bosons in a double-well potential is investigated 
from the point of view of many-body quantum mechanics in the framework of the two-mode model. 
For small initial population imbalances, coherent superpositions of almost equally spaced energy 
eigenstates lead to Josephson oscillations. The suppression of tunneling at population imbalance 
beyond a critical value is related to a high concentration of initial state population in the region of 
the energy spectrum with quasi-degenerate doublets resulting in imbalance oscillations with a very 
small amplitude. For unaccessible long times, however, the system recovers the regime of Josephson 
oscillations. 
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The understanding of many-body quantum systems 
from the theoretical and experimental points of view 
has undergone a considerable development during the 
past decade. Unifying concepts of several branches of 
physics are under development, creating an interdisci- 
plinary scenario for the understanding of quantum me- 
chanical paradigms. One of the simplest many-body sys- 
tems to be realized experimentally and studied theoreti- 
cally are ultracold bosons in a double- well potential. This 
system is very rich exhibiting a great variety of quan- 
tum phenomena such as interference tunneling/self- 
trapping 2, 3, A ^5, 6, 3|i entanglement of macroscopic 
superpositions Lately this system has been exten- 
sively studied, especially after the implementation of sev- 
eral experiments in the area. The usual theoretical ap- 
proach to weakly interacting Bose-Einstein condensates 
(BECs) is the mean-field approximation, a nonlinear 
Gross-Pitaevski equation [llETlO, 11, IzlElElElEl, 
which has proven very adequate in explaining a wide va- 
riety of experiments. 

More recently, the dynamics of population distribu- 
tion between two or more wells of an optical lattice 
have been experimentally investigated. In particular, 
Josephson oscillations have been observed in a ID op- 
tical lattice 17j.2fi] a-nd recently the density distribution 
of tunneling particles is directly observed 0|. In 

this experiment, initial population diff'erences between 
the left and right well components are controlled by 
loading the BEC into an asymmetric double-well poten- 
tial. The Josephson dynamics is initiated at t = by 
non-adiabatically changing the potential to a symmet- 
ric double-well. When the initial population imbalance 
is below a critical value, the system presents Josephson 



oscillations between the two sides of the well. However, 
above this critical value tunneling is not observed. Based 
on a mean field treatment, this is usually attributed to 
macroscopic self-trapping. In the present work, we dis- 
cuss an alternative approach to this system based on ex- 
act numerical solutions of the two-mode Bose-Hubbard 
Hamiltonian [l9j |: 
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where a| and ai are, respectively, the creation and an- 
nihilation operators of a boson in the i-th site, the 
occupation number in the i-th site, U is the two-particle 
on-site interaction, S is the trap depth difference between 
the two wells, J is the tunneling matrix element between 
adjacent sites which characterizes the strength of the 
tunneling term. 

In this Letter we present a different view on the experi- 
mental non-observation of tunneling for initial population 
imbalances larger than a critical value Zc- The dynamics 
of the system is related to the properties of the energy 
spectrum and to the distribution of populations and co- 
herences of the initial wave-function. It is important to 
notice that this model displays many-body features such 
as highly correlated many-body eigenstates. In these 
eigenstates the number of atoms in each well is not de- 
fined but features correlated fluctuations. According to 
the present model, for certain range of the parameter J/U 
the spectrum of the many-body Hamiltonian presents 
a regk)n with nearly degenerate doublets, as shown in 
Ref. p] and Fig.^ For sufficiently large initial population 
imbalance, mainly doublets will participate in the tunnel- 
ing process. The tunneling period will be approximately 
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given by t ~ ?i/ AEdoubiets where Ai?doubict stands for the 
small inner energy difference of the doublets. This energy 
splitting as Ai;doubiet/C/ = [2N {J /U)^{N -1)1 (calcu- 
lated via the doorway method of Ref [23j). For example, 
for N = 100 particles, J/U = 0.333, A£;doubict/C/ 
3.76 X lO^^''^ 0. This energy splitting is usually 
very small and hence will result in periods which are 
much longer than the time window of the experiments. 
In the experimental situation of Ref. only much 
shorter times are observed which is far from the time 
the tunneling actually happens, giving the idea of "self- 
trapping". This short time regime is characterized by 
oscillations with small amplitude which are related to 
the small coherence between the two pairs of adjacent 
quasi-degenerate doublets of the spectrum. For this case, 
the oscillation period is given by the inverse of the en- 
ergy splitting of adjacent quasi-degenerate doublet, e.g., 
t = h/AE, with AE/U ~ (iV - 1). For this reason, the 
interpretation of the non-observation of tunneling in the 
model is radically different from the one coming from the 
Gross-Pitaevski equation, attributed to the non linearity 
of the mean-field approach The same analysis can be 
extended to the improved two-mode model developed in 
Ref. 13]. 

To illustrate our ideas, we study a hundred particles 
in a double-well potential. The upper graphic in Fig. ^ 
shows the energy spectrum for the symmetric case at 
J/U ~ 3.333 and N ~ 100, which corresponds to the 
parameter A = = 15 used in the experimental situ- 
ation of Ref. [3- The spectrum is composed of two re- 
gions: one region of almost equidistant energy levels and 
one region consisting of quasi-degenerate doublets, (see 
also Ref. Q)- The dynamics of population in these two 
regions differs vastly. It has been shown in Ref. that 
the eigenstates of the system confirm the physical pendu- 
lum characteristics of the eigenstates. The ground state 
is a minimum uncertainty wave packet in both number 
and phase which is centered at the origin. The harmonic- 
oscillator-like low-lying excited states are the analog of 
pendulum librations, and the higher- lying cat like states 
are the analog of pendulum rotor motions, with a clear 
signature of the quantum separatrix state where the li- 
bration and rotation states separate. This separatrix di- 
vides the energy spectrum in two regions, with or without 
quasi-degenerate doublets which can be directly observed 
in the upper graph of Fig. ^ The existence of these two 
distinct spectrum regions has a direct consequence for 
the dynamics. 

The tunneling dynamics is characterized by the time 
evolution of the population imbalance 
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FIG. 1: Upper graph: Hamiltonian energy spectrum for 100 
bosons in a symmetric double-well potential for J/U = 3.333 
corresponding to A = 15. Lower graph: Density plot of the 
occupation probability of symmetric Hamiltonian eigenstates 
|m) as a function of the initial population imbalances z{t = 0) 
for the same parameters as the upper graph. The gray scale 
gives the population of the corresponding eigenstate. Below 
a critical value Zc « 0.5, only states with almost equal energy 
spacing are populated, while above Zc the quasi-degenerate 
doublet states are occupied. 



where \ipit)) = Cme"'^"*!™) and z^m' = 

c*raCm' {Era\{ni - n2)\Em')/N, with Em and |m) be- 
ing, respectively, the eigenvalues and eigenstates of the 
hamiltonian in the symmetric case (5 = in Eq. iQ), 
Cm — (("^l'0(O))) and |'0(O)) is the initial many-body 
wavefunction. Following the experimental approach of 
Ref. the initial condition is then prepared by non- 
adiabatically transferring the the ground state of an 
asymmetric double-well potential 6 =/= to the symmet- 
ric double-well. This corresponds to a projection of the 
ground state wavefunction for the asymmetric Hamilto- 
nian (5 7^ to the basis \m). The initial state thus crit- 
ically depends on the value of the parameter S which is 
characterized by the initial population imbalance z(0). 
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Two regions of the energy spectrum can be accessed 
choosing appropriate initial conditions for the many- 
body wavefunction of the system. In the lower graph 
of Fig. n the occupation probability |cm|^ of the eigen- 
states \m) as a function of the initial population imbal- 
ance 2:(0) is depicted. For small z(0) only the lowest few 
energy eigenstates are populated. The spread of the dis- 
tribution of populated states increases as z(0) increases, 
but for ^(0) below a critical value Zc ~ 0.5 only energy 
eigenstates with nearly equal energy splitting are occu- 
pied (see upper graph in figure^]). As follows from Eq.|21 
coherences between these states will lead to a macro- 
scopic oscillation of the population imbalance at a fre- 
quency given by the energy splitting. For larger values 
z{Q) > Zc this scenario changes and only doublet states 
are occupied. These states are characterized by superpo- 
sition states in the basis of states confined to the left and 
right side of the double. As these two states are almost 
degenerate, coherences among them do not contribute to 
the oscillatory term in Eq. |21 

The dynamical behavior of the system z(t) crucially de- 
pends on the competition of the two ingredients of z„ „j, 
the correlation matrix of the initial state in the basis of 
the eigenstates of the symmetric hamiltonian cj^c„ and 
the correlation matrix of the observable, the population 
imbalance in the basis of the eigenstates of the symmet- 
ric hamiltonian {Em\ini — n2)\E„). The second term is 
responsible for the suppression of the population oscilla- 
tions of the first term. As the occupied eigenstates are 
almost equally spaced for z(0) < Zc, the population im- 
balance z(t) oscillates around zero at a fundamental oscil- 
lation frequency given by the mean energy spacing of the 
occupied states, which is given by the plasma frequency 
ujp = 2JVTTA, with A = {UN/2J) p. This behavior 
constitutes the regime of Josephson oscillations as shown 
by the left graph of the upper curves in Fig. |5| Small 
deviations from equal energy spacing lead to a damping 
of the oscillation on a time scale inversely proportional 
to the number of occupied states and the corresponding 
differences in energy splitting. The oscillations undergo 
revivals on time scales given by the inverse of the fre- 
quency difference of adjacent eigenstates. The revivals 
can be seen in the right graph of the upper curves in 
Fig. 13 For larger values of z(0) > Zc{0), where doublet 
states are occupied with negligible coherences between 
adjacent doublet, one observes two different time scales, 
characterized by a regime with small or large oscillation 
amplitude. For the smaller time scales, the population 
imbalance is locked to its initial value with small resid- 
ual oscillations as shown by the left graph of the lower 
curves in Fig. |21 The frequency of these small oscilla- 
tions is determined by the inverse of the energy splitting 
of two quasi-degenerate doublets, i.e. AE ~ NU. The 
latter regime is commonly the " Self- Trapping" regime 
observed in Ref.Q- It should however be noted, that on 
larger time scales T ~ ?i/ Ai^doubict, which are far beyond 



experimental observation, the atoms still undergo col- 
lective tunneling resulting in oscillatory behavior of the 
population imbalance around zero, since the initial condi- 
tion guarantees a high occupation of the quasi-degenerate 
doublets, as one can see in the lower graph of figure ^ 

The time average of the population imbalance z{t) over 
50 plasma periods {2Tr/ujp) is shown in Fig|31for different 
values of A as a function of the initial population imbal- 
ance z(0) and gives the critical value Zc for which the 
system crosses from the Josephson regime (oscillations 
around z = 0) to a self-trapping regime. From Fig|31one 
sees that as the parameter A decreases, the critical value 
Zc increases for a fixed number of particles N. This can 
be understood in terms of the structure of the symmet- 
ric hamiltonian energy spectrum. The number of quasi- 
degenerate doublets increases as the tunneling parame- 
ter J/U decreases (A increases). As a consequence, the 
quasi-degenerate doublets have lower energies which can 
be accessed by lower values of Zc, providing an increase 
of the self-trapping region, since the region of high distri- 
bution of the initial population imbalance in the quasi- 
degenerate doublets also increases. The critical values 
Zc of Fig 121 follow very well the semiclassical prediction 
Zc ~ 2Vl + A/A. 

In spite of the obvious advantages of an exact ap- 
proach, there are several limitations to the model. The 
two mode assumption is an approximation to the real 
situation. More restrictive, however, is the hypothesis 
that the coefficients J and U remain constant in time. 
This approximation is only valid when the many-body 
interactions produce small modifications on the ground 
state properties of the individual wells. This is true if the 
on-site interaction energy is much smaller than the level 
spacing of the external trap (see In this case the 

number of atoms N <C where a is the scattering 

length and tq is the position uncertainty in a harmonic 
oscillat or g round state. In the experiment with of 
ref. 13, lU N ^ 200. However, qualitative behavior is 
still observed. 

We have investigated the dynamics of the imbalance 
population of bosons in a double-well potential from the 
point of view of a many-body Hamiltonian in the frame- 
work of the two-mode model. Although the model is not 
realistic enough for large numbers of particles (TV > 200), 
it points out a completely different explanation for the 
suppression of tunneling. There are no nonlinearities in 
this system and therefore the tunneling process is ex- 
plained solely in terms of initial conditions and spectral 
properties of the many-body Hamiltonian. In this con- 
text, the spectrum is divided in two regions, one of them 
consisting of quasi-degenerate doublets. When the ini- 
tial condition is small such that z(0) < Zc, the occupa- 
tion probability is larger in the lower part of the energy 
spectrum with quasi-equidistant levels, contributing with 
large oscillation amplitude with a plasma frequency ujp. 
When the initial population imbalance is large enough 
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same analysis can be extended to the improved two-mode 
model of ref 0| , which includes the terms reflecting the 
physics of the transfer of atoms from one well to the other 
due to collision excluded from the two-mode model, since 
the energy spectrum of the improved model has the same 
structure as the two-mode model, i.e. it is also divided 
into a region with and without quasi-degenerate doublets. 
Therefore, the same qualitative physical picture would 
result by this extension of the present model. 
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FIG. 2: Time evolution of the population imbalance 100 par- 
ticles for J/U = 3.333 and initial conditions z{Q) = 0.3 (upper 
graph) and z{0) — 0.7 (lower graph). 
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FIG. 3: Time average of z{t)/z{0) over 50 plasma periods 
as a function of the initial population imbalance (^(0)) for 
100 particles. Josephson oscillations result in z{t) — while 
z{t) 7^ indicates the regime of suppressed tunneling. The 
parameter A — NU/2J determines the critical population 
imbalance separating these two regimes. 

z(0) > Zc only the quasi-degenerate doublets are appre- 
ciably occupied. For small time range, the system os- 
cillates with a very small amplitude, whose oscillation 
period is given t ~ fi/{U{N — 1)). However, for longer 
times, the system recovers the Josephson oscillation be- 
havior, with a large amplitude and oscillatory period of 

t ~ ASdoubict/C/ = C^^N-iy.^ ) t^i^ tunnel- 
ing time given for a large number of particles is beyond 
observational possibility. Hence, the "self-trapping" ob- 
served in the experiment of ref.Q is related to the be- 
havior of the imbalance population for short time. The 
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